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Abstract 



In the paper the closed Friedmann-Robertson-Walker model with quantization in the presence of the 
positive cosmological constant and radiation is studied. For analysis of tunneling probability for birth of 
an asymptotically deSitter, inflationary Universe as a function of the radiation energy a new definition of a 
"free" wave propagating inside strong fields is proposed. On such a basis, tunneling boundary condition is 
corrected, penetrability and reflection concerning to the barrier are calculated in fully quantum stationary 
approach. For the first time non-zero interference between the incident and reflected waves has been taken 
into account which turns out to play important role inside cosmological potentials and could be explained 
by non-locality of barriers in quantum mechanics. Inside whole region of energy of radiation the tunneling 
£Nl . probability for the birth of the inflationary Universe is found to be close to its value obtained in semiclassical 

^ ' approach. The reflection from the barrier is determined for the first time (which is essentially differs on 1 at 

the energy of radiation close to the barrier height). The proposed method could be easily generalized on the 
cosmological models with the barriers of arbitrary shape, that has been demonstrated for the FRW-model 
with included Chaplygin gas. Result is stable for variations of the studied barriers, accuracy are found to 
be 11-18 digits for all coefficients and energies below the barrier height. 
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1 Introduction 

In order to understand what really happens in the formation of the Universe, many people came to the point 
of view that a quantum consideration of this process is deeper. The first papers with the quantum approach 
for the description of Universe formation and its initial expansion may be [TJ[3], and shortly afterwards many 
other papers appeared in this field, pointing to a rapid development of the quantum approach in cosmology (for 
example, see the first some papers [3HT0] and some discussions in Refs. [TT][T2] with references therein). 

Today, among all variety of models one can select two approaches which are the most prevailing: these are 
the Feynman formalism of path integrals in multidimensional spacetime, developed by the Cambridge group 
and other researchers, called the "Hartle-Hawking method" (for example, see Ref. [3]), and a method based 
on direct consideration of tunneling in 4-dimensional Euclidian spacetime, called the "Vilenkin method" (for 
example, see Refs. [3][8]|9llll])- Here, according to Ref. [13], in the quantum approach we have the following 
picture of the Universe creation: a closed Universe with a small size is formed from "nothing" (vacuum), where 
by the word "nothing" one refers to a quantum state without classical space and time. A wave function is used 
for a probabilistic description of the creation of the Universe and such a process is connected with transition of 
a wave through an effective barrier. Determination of penetrability of this barrier is a key point in estimation 
of duration of the formation of the Universe, and dynamics of its expansion in the first stage. 

However, in majority of these models (with the exception of some exactly solvable models) tunneling is 
practically studied in details in the semiclassical approximation mainly (for example, see Refs. |11II12| ). An 
attractive side of such an approach is its simplicity in the construction of decreasing and increasing partial 
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solutions for the wave function in the tunneling region, an outgoing wave function in the external region, and 
a possibility to define and to estimate in a simply enough way a penetrability of the barrier, which can be 
used for obtaining the duration of the nucleation of the Universe. A tunneling boundary condition |11[|13) could 
seems to be the most natural and clear, where the wave function should represent an outgoing wave only in 
the enough large value of the scale factor a. However, whether is such a wave really free in the asymptotic 
region? If to draw attention on increasing of a modulus of the potential with increasing of the scale factor a 
and increasing of a gradient of such a potential, used with opposite sign and having a sense of force, acting 
"through the barrier" on this wave, then one come to a serious contradiction: influence of the potential on 
this wave increases strongly at increasing of a\ Now a new question has been appeared: what should the wave 
represent in a general in the cosmological problem? This problem connects with another (new) else more general 
one in general quantum physics — a real importance to define a "free " wave inside strong fields and we need 
in mathematical stable tools for its study and work. It becomes unclear whether a connection between exact 
solutions for the wave function at turning point and "free" wave defined in the asymptotic region is corrected. 
If in frameworks of semiclassical approach (up to the second order) it is corrected then any deformations of tail 
of the potential after the barrier cannot change the penetrability absolutely, while force mentioned above and 
acting on the asymptotic "free" wave could be increased up to infinity! Answer on such misunderstanding could 
be found in non-locality of definition of the penetrability in quantum mechanics, which is reduced to minimum 
in the semiclassical approach (i. e. this is so called "error" of the cosmological semiclassical approach). Now, if 
to refuse application of the semiclassical approximation, then we come to a necessity to define such a wave as 
maximal correctly and accurately as a possible. The problem of the correct definition of the wave in cosmology is 
reinforced else more, if to calculate incident and reflected waves in the internal region. Even with the known exact 
solution for the wave function there is uncertainty in determination of these waves! But, namely, the standard 
definition of the coefficients of penetrability and reflection is based on them. In particular, I have not found 
papers yet where the coefficient of reflection is defined and estimated in this problem (which differs essentially 
from unity at the energy of radiation close to height of the barrier and, therefore, such a characteristics could 
be interesting from a physical point of view). Note that the semiclassical approximation prejudices a possibility 
of its definition at all [14] . 

Thus, in order to estimate probability of the formation of the Universe as accurately as possible, we need 
in the fully quantum definition of the wave. Note that the non-semiclassical penetrability of the barrier in the 
cosmological problems has not been studied in details and, therefore, a development of fully quantum methods 
for its estimation is a perspective task. I would like to note undoubtedly perspective direction of researches [151 , 
where the penetrability was estimated on the basis of tunneling of wave packet through the barrier. However, 
a stationary boundary condition has uncertainty that could lead to different results in calculations of the 
penetrability. The stationary approach could allow to clarify this deeply. It is able to give stale solutions for 
the wave function (and results in Ref. [16] have confirmed this at zero energy of radiation), uses the standard 
definition of the coefficients of the penetrability and reflection, is more accurate in their estimation. 

Aims of this paper are: (1) to define the wave in the quantum cosmological problem; (2) to construct the 
fully quantum (non-semiclassical) stationary method of determination of the coefficients of penetrability of the 
barriers and reflection from them on the basis of such a definition of the wave; (3) to estimate how much 
the semiclassical approach is differed in the estimation of the penetrability from the fully quantum one. In 
order to achieve this, at first it needs to construct tools for stable calculation of partial solutions of the wave 
function. In order to resolve the questions pointed out above, we shall restrict ourselves by a simple cosmological 
model, where the potential has a barrier and internal above-barrier region. The paper is organized so. After 
a short description of closed Fricdmann-Robertson- Walker (FRW) model with quantization in a presence of 
a positive cosmological constant and radiation in Sec. 2, in Sec. 3 a new formalism for calculation of two 
linear independent partial solutions for the wave function of the Universe for the scale factor inside the region 
< a < 100 and the energy of radiation from zero up to the barrier height is presented. In Sec. 4 a fully quantum 
definition of the wave is formulated (for the first time), a new approach for determination of the incident, 
reflected and transmitted waves relatively the barrier is constructed on such basis, the boundary condition is 
corrected. In Sec. 5 a fully quantum stationary method of determination of coefficients of penetrability and 
reflection relatively the barrier with analysis of uniqueness of solution is presented, where at first time non-zero 
interference between the incident and reflected waves is analyzed and for its estimation the coefficient of mixing 
is introduced. In such an approach the penetrability and reflection for the barrier for the studied cosmological 
model with A — 36, B ~ 12 A parameters at A = 0.01 are calculated with estimation of accuracy and in 
comparison with the semiclassical results and with results of non-stationary quantum approach [15] . In Sec. 6 
the penetrability is estimated for the FRW-model with included Chaplygin gas which has internal hole before 
the barrier. Conclusions finalize the paper. 
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2 Cosmological model in the Friedmann— Robertson— Walker metric 
with radiation 



Let us start from a case of a closed (k = 1) FRW model in the presence of a positive cosmological constant 
A > and radiation. The minisuperspace Lagrangian has the form (see Appendix A, also Ref. |13) . (11), p. 4): 

£(a ' d) = 8^(- fl2+ ^^ a2p(a) )' p{a) = PA + ^wy pa = s7g' (1) 

where a is scale factor, a is derivative a with respect to time coordinate t, p (a) is a general expression for 
the energy density, p ra d(d) is component describing the radiation in the initial stage (equation of state for 
radiation is p(a) = p r ad(d)/3, p is pressure). The passage to the quantum description of the evolution of the 
Universe is obtained by the standard procedure of canonical quantization in the Dirac formalism for systems 
with constraints. In result, we obtain the Wheeler-De Witt (WDW) equation (see Ref. [13], also [Tl[2l[22]). 
which after multiplication on factor and passage of the item at the component with radiation p ra( j into right 
part transforms into the following form (see Appendix A): 

{-£ + v(.)}*.)-w.), n«) = (i^)%<. a -f^, < 2 » 

where <p(a) is wave function of Universe. This equation looks similar to the one-dimensional stationary 
Schrodinger equation on semiaxis (of the variable a) at energy E^-ad with potential V (a). For further anal- 
ysis it is convenient to use the system of units where 8ttG = M p = 1, and to rewrite V (a) in a generalized 
form: 

V(a) = Aa 2 - Ba 4 . (3) 

In particular, for the Universe of the closed type (k = 1) we obtain A = 36, B = 12 A (this potential coincides 
with [15]). 

In order to find wave function we shall needs to know shape of the potential close to turning points. Let us 
find the turning points dt p ,i n an d fltp, out concerning the potential ([3]) at energy £" ra d: 



dtp, in — 




Let us expand the potential V(a) ^ in powers g ou t = a — a tp (where as at p the point a tPi i n or a tp . ut is used, 
close to which we find expansion), where (for small q) we restrict ourselves by the linear item only: 

V(q) = V + V 1 q, (5) 

where the coefficients Vq and V\ are: 



V = V(a = dtp, m) = V(a = a tp , out) = Aa 2 t ' - B a 4 = E r 



rad ; 



T/ (out) _ , A f ABE r&d \ / / ABE iad \ 

^ = - 2A ^2By-^)\ l + V-^y (6) 



Vr - 2 ..^(l-^)(l-^l-^ 

Now eq. ([2]) transforms into a new form at variable q with potential V(q): 

d 2 



dq 



2 ip(q) + Vi q <p(q) = 0. (7) 



3 Wave function: its behavior and partial solutions 

The wave function is known to oscillate above the barrier and increase (or decrease) under the barrier without 
any oscillations. So, in order to provide a linear independence between two partial solutions for the wave 



3 



function effectively, I look for the first partial solution increasing in the region of tunneling and the second 
one decreasing in this tunneling region. At first, I define each partial solution and its derivative at a selected 
starting point, and then I calculate them in the region close enough to this point using the method of beginning 
of the solution presented in Appendix IB. II Here, for the partial solution which increases in the barrier region, 
as the starting point I use internal turning point et tPi j n at non-zero energy -E ra d or equals to zero a — at 
null energy -E ra d, and for the second partial solution which decreases in the barrier region, I select the starting 
point to be equal to external turning point a tp . ut- Then both partial solutions and their derivatives I calculate 
independently in the whole required range of a using the method of continuation of the solution presented in 
Appendix IB. 21 which is improvement of the Numerov method with a constant step. By such a way, I obtain 
two partial solutions for the wave function and their derivatives in the whole studied region. 

In order to clarify how much the proposed approach gives convergent (stable) solutions, for a comparison 
we shall use paper [15] with the published wave function (see (9) p. 5, I use the potential in eq. <j3j) ) - Let us 
consider a behavior of the wave function. The first partial solution for the wave function and its derivative 
in my calculation are presented in Fig. [TJ which increase in the tunneling region and have been obtained at 
different values of the energy of radiation E^ad . Form these figures one can see that the wave function satisfies 




Figure 1: The first partial solution for the wave function and its derivative at different values of the energy of radiation 
_B ra d, increasing in the tunneling region (curve 1, blue, is for the wave function; curve 2, green, for the derivative of this 
wave function): (a) £ rad = 10; (b) £ rad = 1000; (c) £ rad = 2000 

the rules of behavior of the wave functions inside sub-barrier and above-barrier regions |24j . Starting from 
very small a, the wave function has oscillations and its maximums increase monotonously with increasing of 
a, which corresponds to behavior of the wave function in the internal region before the barrier (this becomes 
more obvious after essential increasing of scale, see left panel in Fig. [J). Further, with increasing of a the wave 
function increases monotonously without any oscillations, that points out on transition into the tunneling region 
(one can see this in a logarithmic presentation of the wave function, see central panel in Fig. [5]). A boundary 
of such a transformation in behavior of the wave function must be the point of penetration of the wave into 
the barrier, i. e. the internal turning point dtp, in- Further, with increasing of a the oscillations are appeared in 
the wave function, which could be possible inside the above barrier region only (in the right panel of Fig. [5] one 
can see that such a transition is extremely smooth that characterizes the accuracy of the method positively) . 
A boundary of such a new transformation in the behavior of the wave function must be the point of leaving of 
the wave from the barrier outside a tp , ou t ■ Like Ref. |16] , but at arbitrary non-zero energy -E ra d in the external 
region, starting from Ot p , out , I obtain again a monotonous increasing of maximums of the derivative of the wave 
function and a smooth decreasing of this wave function. One can see that the derivative is larger essentially 
than the wave function. At enough large values of a, i. e. in the region, which can be called asymptotic one, I 
obtain the smooth continuous solutions, achieving a = 100 (in Ref. [T5] the maximal value is a = 30). 

The second partial solution of the wave function and its derivative in my calculation at different values of 
the energy of radiation -E ra d is presented in next Fig. [3] which decrease in the region of tunneling. According to 
analysis, this solution close to the turning points, in the tunneling region, in the sub-barrier and above-barrier 
regions looks like the first partial solution, but with a difference that now maximums of the wave function and its 
derivative are larger essentially in the external region in a comparison with the internal region, and amplitudes 
in the tunneling region are decreased monotonously. 

Comparing all these pictures of the wave function with the results of Ref. [15], one can see that the wave 
function in such an approach is essentially more continuous, has no any divergencies and its behavior is clear 
anywhere. From here I conclude that in determination of the wave function and its derivative at arbitrary 
energy of radiation the developed method is essentially more quick, stable and accurate in a comparison with the 
non- stationary quantum approach in Ref. fl5f . Note the following. 
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Figure 2: The first partial solution for the wave function and its derivative at the energy of radiation -E ra d = 2000 (curve 
1, blue, is for the wave function; curve 2, green, for the derivative of this wave function) 
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Figure 3: The second partial solution for the wave function and its derivative at different values of the energy of radiation 
Srad, decreasing in the tunneling region (curve 1, blue, is for the wave function; curve 2, green, for the derivative of this 
wave function): (a) E lad = 10; (b) S rad = 1000; (c) £ rad = 2000 



• With increasing a, a period of the oscillations both for the wave function and its derivative decreases 
uniformly in the external region and increases uniformly in the internal region (it partially obtained 
earlier in Ref. [16 at i? ra d = 0). 

• At larger distance from the barrier (i. e. at increasing of a in the external region and at decreasing of a in 
the internal region) it becomes more difficult to achieve the convergent continuous solutions for the wave 
function and its derivative (it partially obtained earlier in Ref. [16 at -E ra d = 0). 

• A number of oscillations of the wave function in the internal region increases with increasing of the energy 
of radiation E la d (a new opened property) . 



4 A fully quantum definition of "free" wave inside strong fields 

4.1 A problem of definition of "free" wave in cosmology and correction of the 
boundary condition 

Which boundary condition should be used for such a wave function, which describes leaving of the wave from the 
barrier outside the most correctly and accurately? A little variation of the boundary condition leads to change 
of the fluxes concerning the barrier and, as result, it changes the coefficients of penetrability and reflection. So, 
a proper choice of the boundary condition is extremely important. However before, let us analyze how much 
the choice of the boundary condition is natural in the asymptotic region. 

• In tasks of decay in nuclear and atomic physics such potentials of interactions are used, which tend 
to zero in the asymptotic region. In these tasks an application of the boundary condition at limit of 
infinity does not give questions. In cosmology we deal with another, principally different type of the 
potential: with increasing of the scale factor a modulus of this potential increases. A gradient of the 
potential, used with opposite sign and having a sense of force acting on the wave, increases also. Therefore, 
here there is nothing mutual with free propagation of the wave in the asymptotic region. Thus, a direct 
passage of the application of the boundary condition in the asymptotic region into cosmological problems 
looks very questionable. 
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• Results in Ref. [16 reinforce a seriousness of this problem, which show that with increasing of the scale 
factor a the region is enlarged where solutions for the wave function (and its two partial solutions) are 
stable. According to Ref. [IB], the scale factor a in the external region is larger, the a period of oscillations 
of each partial solution for the wave function is smaller . This requires with increasing of a to decrease 
continuously step of calculations and to increase time of calculations of the wave function for keeping the 
same accuracy. This increases errors in calculation of the wave (if to set it to be free in the asymptotic 
region), and in result the final found solution can have nothing mutual with proper one! From here a 
natural conclusion follows on a principal impossibility to use practically the boundary condition at infinity 
for calculation of the wave (in supposition if we know it maximally accurately in the asymptotic region), 
if we like to pass from the semiclassical approach to the fully quantum one. Principally another situation 
exists in problems of decay in nuclear and atomic physics where calculations of the wave in the asymptotic 
region are the most stable and accurate! 

• One can add a fact that it has not been known yet whether the Universe expands at extremely large value 
of the scale factor a. Just the contrary, it would like to clarify this from a solution of the problem, however 
imposing a condition that the Universe expands in the initial stage. 

On such a basis, I shall introduce the following definition of the boundary condition: 

The boundary condition should fix the wave function so that it represents the wave, interaction 
between which and the potential barrier is minimal at such a value of the scale factor a where action 
of this potential is minimal as possible. 

A propagation of the wave defined by such a way will be closed maximally to free one for the used potential 
and at used value of the scale factor a (I call such a wave conditionally "free"). However, with a purpose to 
give a mathematical formulation for this definition we are confronted with two questions, which we will have to 
resolve: 

1. What should the free wave represent in general in a field of the cosmological potential of arbitrary shape? 
How could it be defined the most correctly in enough small neighborhood of arbitrary selected coordinate? 

2. At which coordinate is the application of such a boundary condition the most corrected? 

At first, let us consider the second question where it needs to impose the boundary condition on the wave 
function. Which should this point be: or this is a turning point (where the potential coincides with energy of 
radiation), or this is a coordinate where a gradient from the potential (having a sense of force of interaction) 
becomes zero, or this is a coordinate where the potential becomes zero? With a purpose to define the wave 
with free propagation, we are looking for such a coordinate where the action by the potential on this wave is 
minimal. 

We define this coordinate where the force acting on the wave is minimal. Here, we define the force 
as the gradient of the potential used with opposite sign. 

It turns out that according to such a (local) definition the force in the external region is minimal at the external 
turning point at Pj0 ut- Also, the force, minimally acting on the wave incident on the barrier in the internal region 
and on the wave reflected from it, is minimal at the internal turning point <2tp,in- Thus, we have just obtained 
the internal and external turning points where we should impose the boundary conditions in determination of 
the waves. At zero energy null of the potential coincides with the turning point and therefore a choice of the 
point of leaving of the wave from the barrier outside in imposing of the boundary condition in Ref. |16j was 
made properly. Here, it could be interesting to analyze a condition of equality to zero of the gradient of the 
potential. But this point is located in the region of the barrier where there is tunneling and, therefore, we shall 
not study this case in this paper. 

4.2 Definition of the wave minimally interacting with the potential 

Now we shall be looking for a form of the wave function in the external region, which describes maximally 
accurately the wave, propagation of which is the closest maximally to "free" one in the external region at the 
turning point Ot p , out and is directed outside. Let us return back to eq. ([7]) where the variable q — a — at Pj out 
has been introduced. Changing this variable to 

e= | Fi (out)|l/3^ (g) 



<» 



this equation is transformed into 

|U(0+f¥>(0 = 0. ( 9 ) 

From quantum mechanics we know two linearly independent exact solutions for the function ip(£) in this 
equation — these are the Airy functions Ai (£) and Bi (£) (for example, see Ref. [25], P- 264-272, 291-294). 
Expansions of these functions into power series at small £, their asymptotic expansions at large |£|, their 
representations through Bessel functions, zeroes and their asymptotic expansions are known. We have some 
integrals of these functions, and also the form of the Airy functions in the semiclassical approximation (which 
can be applied at large |£|). In some problems of the analysis of finite solutions <^(£) in the whole range of £ 
it is convenient to use the integral representations of the Airy functions (see eq. (10.4.32) in Ref. [25], p. 265, 
a = 1/3, taking into account sign in eq. (10.4.1)): 

// :! 







1 


/ 


7T 











+ OC 


1 


/ 


7T 










Ai (±£) = - / cos ( y T & ) du, 



Bi(±0 = - 



(10) 



ex P I — g- T J + sin I — =F 



du. 



Furthermore, we shall be interested in the solution of the function which describes the most accurately the 
outgoing wave in the range of a close to the at p point. However, it has been not clear what the wave represents 
in general near the point at p in the potential studied, and which linear combination of the Ai(£) and Bi(£) 
functions defines it in the most accurate way. 

The clearest and most natural understanding of the outgoing wave is given by the semiclassical consideration 
of the tunneling process. However, at the given potential the semiclassical approach allows us to define the 
outgoing wave in the asymptotic region only (while we can join solutions in the proximity of at p by the Airy 
functions). But it is not clear whether the wave, defined in the asymptotic region, remains outgoing near the 
atp. During the whole path of its propagation outside the barrier the wave interacts with the potential, which 
must inevitably lead to a deformation of its shape (like to appearance of a phase shift in the scattering of a 
wave by a radial potential caused by interaction in scattering theory). Whether does it turn out the potentials 
used in cosmological models give a significantly larger change of the shape of the wave caused by interaction in 
a comparison with the potentials used, for example, for the description of nuclear collisions in the framework 
of scattering theory? Moreover, for the given potential there is a problem with obtaining convergence in the 
calculation of the partial solutions for the wave function in the asymptotic region. According to our calculations, 
a small change in the range of the definition of the wave in the asymptotic region leads to a significant increase 
of errors, which requires one to increase the accuracy of the calculations. Therefore, we shall be looking for a 
way of defining the outgoing wave not in the asymptotic region, but in the closest vicinity of the point of escape, 
a tp . In a search of solutions close to the point a tp , i. e. at small enough |£|, the validity of the semiclassical 
method breaks down as |£| approaches zero. Therefore, we shall not use the semiclassical approach in this paper. 

Assuming the potential V(a) to have an arbitrary form, we define the wave at the point at p in the following 
way. 

Definition 1 (strict definition of the wave) The wave is such a linear combination of two partial solutions 
of the wave function that the change of the modulus p of this wave function is closest to constant under variation 
of a: 

d 2 



da* Pia) 



0. (11) 



According to this definition, the real and imaginary parts of the total wave function have the mutually closest 
behaviors under the same variation of a, and the difference between possible maximums and minimums of the 
modulus of the total wave function is the smallest. For some types of potentials (in particular, for a rectangular 
barrier) it is more convenient to define the wave less strongly. 

Definition 2 (weak definition of wave): 



The wave is such a linear combination of two partial solutions of wave function that the modulus p 
changes minimally under variation of a: 

d 



da P{a) 



->0. (12) 

=a tp 
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According to this definition, the change of the wave function caused by variation of a is characterized mainly 
by its phase (which can characterize the interaction between the wave and the potential). 

Subject to this requirement, we shall look for the solution for the function <£>(£) in the following form: 

¥>(0=T.*<+>(0, (13) 

where 

Umax 

r 3 

*(±)(0 = / exp±i(-y + /(Ou)du. (14) 
o 

where T is an unknown normalization factor, /(£) is an unknown continuous function satisfying /(£) — > const 
at £ — > 0, and it ma x is the unknown upper limit of integration. In such a solution, the real part of the function 
/(£) gives a contribution to the phase of the integrand function only while the imaginary part of /(£) deforms 
its modulus. 

Let us find the first and second derivatives of the function 4 r (^) (a prime denotes a derivative with respect 
to £): 

d ¥ ± 'H0 = ±t J /'ti exp±t(-y+/(0«) d«, 



(15) 



d£ 



ii max 

^* (±) (0 = / (±i/"u-(/') 2 « a ) exp±i(-y + /(£)«) du. 
o 

From this we obtain: 

u max 

^* (±) (0+e* (±) (0= / (±i/"t*-(/')V+f) exp±i(-y +/(0«) d«. (16) 



Considering the solutions at small enough values of |£|, we represent /(£) in the form of a power series: 

+ OG 
71=0 

where /„ are constant coefficients. The first and second derivatives of /(£) are 

= 

d£ 



"too " = ° (18) 



f'(0 = 4 /(0 = £(n + !) (" + 2 ) € 

^ 71=0 



Substituting these solutions into eq. (|15|) . we obtain 

*(±) (£) + e = j J (± 2m / a - u 2 ffj + (± 6zw / 3 - 4 M 2 A/ 2 + l) £+ 



d£ 2 



f OO r 71 -| j , u 3 . 

±iu(n+ l)(n + 2) / n+2 - u 2 (n - m + l)(m + 1) f n -m+ifm+i f l >exp±i( — 
i=2 L m=o J J v o / 

Considering this expression at small |£|, we use the following approximation: 



du. 

(19) 



w'> 



exp±^ — - + fuj -> exp±i(^ — - + / uJ. (20) 



Then from eq. © we obtain the following condition for the unknown /„ 



: : 2/// /■> - u 2 ffj exp±i (~y + /•>») (/ " + 

+ £ • / (± Gin h ~ 4« 2 hh + l) exp ±» (-y + /o«) d« + 
o 

^max 

+°° /" r 1 / u 3 \ 

+ / ± iu (n + l)(n + 2) /„ +2 - u 2 VVn-m + l)(m + 1) /„- m +i/m+i cxp±i (-— + / u) du = 0. 

(21) 

Requiring this condition to be satisfied for different £ with different powers n, we obtain the following system: 

,3 



J (± 2iu j 2 - u 2 A 2 ) exp ±» (- y + /o«) du = 0, 
o 

Umax 

/ (± 6zu h ~ 4u 2 /i/ 2 + l) exp±i (- y + /o«) du = 0, 
o 

n 3 

±iu (n + l)(n + 2) /„ +2 - u 2 (n - m + l)(m + 1) /„_ m+ i/ m+ i exp±i (-y + /o") du = 0. 



m=0 



(22) 

Assuming the coefficients fo and fi to be given, we find the following solutions for the unknown f$ and 



fn ■ 



f: 



f2 J ±) 
(±) _ , A J 2 f (±) _ , 

3= — ' ~ 7T7> / 3 — 3= 



4f ,(±) A±) j(±) 

4/1/2 J 2 ~ J 



(n - m + l)(m + 1) /m^+i 7 (±) 



P (±) 



2* j{ 



6i Jf 



(±) 



f (±) _ m=0 
Jn+2 ~ 



J 



i (n + l)(n + 2) 



(±) : 



(23) 



where the following notations for the integrals have been introduced: 



1i max ^ Umax ^ lima: 

J^ ±) = J exp±i (-y + / u) du, = J u exp±i (-y + / u) du, J 2 (±) = / 



du. 



(24) 

Thus, we see that the solution (|13l) with taking into account eq. (1T4]) for the function ip (£) has arbitrariness 
in a choice of the unknown coefficients /o, /1 and the upper limit of integration u max . However, the solutions 
found, eqs. (|23|) . define the function /(£) so as to ensure that the equality (fl3|) is exactly satisfied in the region 
of a close to the escape point atp. This proves that the function ip (£) in i/ie /orm, H3\) with taking into account 
eq. \1J$ at an arbitrary choice of fo , f\ and it max is the exact solution of the Schrddinger equation near the 
escape point at p . In order to bring the solution ^(t;) into the well-known Airy functions, Ai (£) and Bi(£), we 
select 

/o = 0, h = l. (25) 

At such a choice of the coefficients fo and /1, the integrand function in the solution up to £ 2 has a constant 
modulus and a varying phase (the coefficient f% deforms the modulus, but it is fulfilled at £ 3 ). Therefore, one 
can expect that the solution ([T^]) at the turning point a tp describes the wave with a maximally proper shape. 

4.3 Total wave function 

Having obtained two linearly independent partial solutions <pi{a) and ^2(0), we make up a general solution (a 
prime is for the derivative with respect to a): 



tp (a) — T ■ (Ci tpi (a) + C 2 <p 2 (a)) , 



Ci = 



iy 2 - *V2 



» ^2 = 
a=a tp , out <Pm-<Pl<P* 



(26) 
(27) 



a— atp, out 
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where T is normalization factor, C\ and Ci are complex constants found from the boundary condition introduced 
above: the <p(a) function should represent an outgoing wave at turning point at P)0 ut- The total wave function 
calculated by such a way for the potential ^ with parameters A = 36, B — 12 A at A = 0.01 at different values 
of the energy of radiation -E ra d is shown in Fig. |4j One can see that a number of oscillations of the wave function 




0,5 1,0 1,5 2,0 

Parameter a 




E =2000 

rad 



Parameter a 




2 3 4 5 6 7 

Parameter a 



Figure 4: The wave function at selected values of the energy of radiation E la _d (curve 1, blue, is for real part of the wave 
function; curve 2, green, for imaginary part of the wave function): (a) -E ra< j = 10; (b) -B ra d = 1000; (c) -Brad = 2000 



in the internal region increases with increasing of the energy of radiation. Another interesting property is the 
larger maximums of the wave function in the internal region at the smaller distances to the barrier for arbitrary 
energy (it has been found for the first time). 

In the next Fig. [5] it has shown how a modulus of this wave function is changed at selected values of the 
energy of radiation. From these figures it becomes clear why the coefficient of penetrability of the barrier is 
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Figure 5: A behavior of the modulus of the wave function at the selected energies of radiation -E ra d: (a) -EWd = 10; (b) 
£ rad = 1000; (c) £ rad = 2000. 

extremely small (up to the energy E xa d — 2000). In order to estimate, how much the boundary condition 
introduced above is effective in construction of the wave on the basis of the total wave function close to the 
external turning point a tp . ut, it is useful to see how the modulus of this wave function is changed close to this 
point. In Fig. [5] the modulus of the found wave function close to the turning points at the energy of radiation 
E ra d = 2000 is shown. Here, one can see that the modulus at at Pl out is constant practically (see left panel in 
Fig. IS]). It is interesting to note that the modulus of so defined wave function close to the internal turning point 
Gtp,in is not changed practically also and is close to maximum (see right panel in Fig. [5]). 

5 A penetrability and reflection in the fully quantum approach 

Let us analyze whether a known wave function in the whole region of its definition allows us to determine 
uniquely the coefficients of penetrability and reflection. 

5.1 Problem of interference between the incident and reflected waves 

Rewriting the wave function <^totai in the internal region through a summation of incident ip- mc wave and reflected 
ip Te { wave: 

Vtotal = Vine + frei, (28) 
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Figure 6: A behavior of the modulus of the total wave function at the energy of radiation £" ra d = 2000 close to the 
turning points (I obtain: atp, in = 8.58, dtp, out = 15.04, also see Table 1): (a) the modulus, decreasing monotonously 
in the tunneling region, with further increasing of a obtains maximums and holes connected with the oscillations of the 
wave function in the external region, but the modulus is not equal to zero (that points out existence of non-zero flux); 
(b) with increasing of a, close to the external turning point a tp ,out the modulus is changed minimally (this demonstrates 
practical fulfillment of the definition for the wave at such a point); (c) transition back close to a tp .m is shown, where 
at increasing of a the modulus with maximums and holes is transformed rapidly into monotonously decreasing without 
maximums and holes, that is connected with transition to the region of tunneling 



-h. 



Jinc H~ Jref H - Jmixed 3 



(29) 



we consider the total flux: 

jOtotal) = i (*>inc + Pref)VU>kc + Pref 

where 

jinc = i (pinc'V'Piac ~ h. C.) , j re { = % (^p re {V(p* ei ~ h. C.) , jmixed = i Cpinc'V <fZet + ¥>ref V</J* nc - h. C.) . 

(30) 

The jmixed component describes interference between the incident and reflected waves in the internal region (let 
us call it as mixed component of the total flux or simply flux of mixing). From constancy of the total flux jtotai 
we find flux j tr for the wave transmitted through the barrier, and: 



jinc = jtr - jref ~ jmixed, jtr = jtotal = COnSt. 



(31) 



Now one can see that the mixed flux introduces ambiguity in determination of the penetrability and reflection 
for the same known wave function. 

5.2 Locality of penetrability, reflection and mixing in the radial task 

In quantum mechanics the coefficients of penetrability and reflection are defined concerning a potential as whole, 
including asymptotic regions. However, in the radial problem of quantum decay such a consideration looks to 
be conditional as the incident and reflected waves should be defined inside internal region from the left of the 
barrier. But, whether does choice of such a region where we use these waves, take influence on the penetrability 
and reflection? In order to formulate these coefficients, we shall include into definitions coordinates where the 
fluxes are defined (denote them as x\ e ft and x r ight): 



T(a;ieft, bright) 



jtr (bright) 
jinc(zieft) 



-R(xieft) 



jrcf(^lcft) 
jinc(^left) 



M(xieft) = 



Jr: 



ed(xieft) 



jinc(zieft) 



(32) 



So, the T and R coefficients determine a probability of transmission (or tunneling) and reflection of the wave 
relatively the region of the potential with arbitrary selected boundaries xi e ft, bright- Tending x r i g ht to the 
asymptotic limit, so defined coefficients should transform into standard ones. From eqs. ([31]) and (|32|) we obtain 
(jtr and j re f are directed in opposite directions, jj nc and jt r — in the same directions): 



\T\ 



\R\-M = 1. 



(33) 



Now we see that condition \T\ + \R\ = 1 has sense in quantum mechanics only if there are no any interference 
between incident and reflected waves which we calculate, and it is enough: 



Jn 



0. 



(34) 
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But a new question has been appeared: whether does this condition allow to separate the total wave function 
into the incident and reflected components uniquely? It turns out that the choice of the incident and reflected 
waves takes essential influence on the barrier penetrability, and different forms of the incident Lp- mc and reflected 
(p rc f waves can give zero flux j m - lx . If we liked to pass from rectangular internal well or some simple presentations 
of waves in semiclassical approach to fully quantum realistic consideration, the problem would become open 
and difficult. Such a situation is typical in quantum cosmology. Whole importance of accurate definition of the 
wave in the quantum cosmological problem becomes clear from here for construction of the total wave function 
on the basis of its two partial solutions, and for the separation of the known wave function into the incident and 
reflected waves in the internal region from the left of the barrier. 

5.3 Wave incident on the barrier and wave reflected from it in the internal region 

One can define the incident wave to be proportional to the function and the reflected wave to be propor- 
tional to the function 

Vtotal(a) = Vine (a) + Vref (a), Vine (a) = / ' * (+) (a) , Vrcf (a) = R ■ * (_) (a) , (35) 

where / and R are new constants found from continuity condition of the total wave function Vtotai and its 
derivative at the internal turning point at p .i n t: 

fr(+) $(-),' _ 

On the basis of these solutions we obtain at the internal turning point atp, int the flux incident on the barrier, 
the flux reflected from it and the flux of mixing. The flux transmitted through the barrier we calculate at the 
external turning point atp,ext- 

5.4 Boundary condition at a = 0: stationary approach versus non-stationary one 

A choice of the proper boundary condition imposed on the wave function is directly connected with questions: 
could the wave function be defined at a = 0, and which value should it be equal to at this point in such a case? 

The complex total wave function is constructed on the basis of its two partial solutions which are calculated 
using methods in Appendix B. Such two partial solutions should be linearly independent, and possible "measure" 
of such linear independence (i. e. nonlinear difference between these two solutions) could define accuracy of the 
final obtained penetrability. In particular, these two partial solutions can be real (not complex), without any 
decrease of accuracy in determination of the total wave function. For any desirable boundary condition imposed 
on the total wave function, such methods should be working. In order to ensure the highest linear independence 
between two partial solutions, I select one solution to be increasing in the region of tunneling and another one 
to be decreasing in this tunneling region. For the increasing partial solution as the starting point a x I use 
internal turning point dtp, in at non-zero energy i? rac i or equals to zero a x — at null energy .E ra d, and for the 
second decreasing partial solution I select the starting point a x to be equal to external turning point dtp, out- 
Such a choice of starting points turns out to give us higher accuracy in calculations of the total wave function 
then if to start calculations of both partial solutions from zero or from only one turning point. So, eq. (173")) and 
furthers in Appendix IB . 1 1 are defined relatively arbitrary non-zero point a x in a general case. These equations 
in Appendixes IB. II and IB. 21 do not absolutely contradict to any desirable boundary condition, and they could 
be extended on the case a x = 0. In latter case the methods in Appendix B should be working at any finite value 
of the total wave function at a = 0. 

In order to obtain the total wave function, we need to connect two partial solutions using one boundary 
condition, which should be obtained from physical motivations. In the light of logic in Introduction, Sect. 4.1 
and 4.2 reinforced by possible interference between the incident and reflected waves in Sect. 5 (which is invisible 
in the semiclassical approach of the second order) it is more natural not to define the total wave function at 
zero (or at infinity) , but to find outgoing wave at such finite value of a in the external region where this wave 
corresponds to our Universe at present time (as this could be connected with possible astronomical observations). 
But, practically, it turns out to define more effectively such a wave at point where forces acting on it are as 
minimal as possible. This is an initial condition imposed on the outgoing wave in the external regiorQ. 

1 For example, using such idea for alpha-decay problem in nuclear physics, we obtain the asymptotic region where the outgoing 
wave is plane (spherical) wave. 



R 



Vtotal 



$(+) 



Vtotal 



(36) 
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Now let us analyze a question: which value could the wave function be equal to at zero? In the paper the 
following basis is used: 

• the wave function should be continuous in the whole spatial region of its definition, 

• we have outgoing non-zero flux in the asymptotic region defined on the basis of the total wave function, 

• we consider a case when this flux is constant. 

The non-zero outgoing flux defined at arbitrary point requires the wave function to be complex and non-zero. 
The condition of continuity of this flux in the whole region of definition of the wave function requires this wave 
function to be complex and non-zero in whole this region. If we included point a = into the studied region, then 
we should obtain the non-zero and complex wave function at such a point also. If we use the basis above, then 
we cannot obtain zero wave function at a = 0. One can pass to nuclear physics where study of such a question 
and its possible resolutions have longer history then in quantum cosmology. As the simplest demonstration, 
one can consider elastic scattering of particles on nucleus (where we have zero radial wave function at r = 0, 
and we have no any divergences), and alpha decay of nucleus (where we cannot obtain zero wave function at 
r = 0). A possible divergence of the radial wave function at zero in nuclear decay problem could be explained 
by existence of source at such a point which creates the outgoing flux in the asymptotic region ( and is a reason 
of such flux). Now the picture becomes clearer: any quantum decay could be connected with source at zero. 
That is why vanishing of the total wave function at a = after introduction of the wall at this point (like in 
Ref. [15]) is not obvious and is only part of all possibilities. 

If we wanted to study physics at zero a — 0, we should come to two cases: 

• If we include zero point into the region of our consideration, we shall come to quantum mechanics with 
included sources. In such a case, condition of constant flux is broken. But more general integral formula 
of non-stationary dependence of the fluxes on probability can include possible sources and put into frame- 
works of the standard quantum mechanics also (see eq. (19.5) in Ref. [2], p. 80). Perhaps, black hole 
could be another example of quantum mechanics with sources and sinks. 

• We can consider only quantum mechanics with constant fluxes and without sources. Then we should 
reduce zero point a = from the region of our consideration. But the formalism proposed in this paper 
remains to be working and is able to calculate the penetrability and reflection coefficients without reduce 
of any accuracy also. 

This could be stationary picture of interdependence between the wave function at zero and the outgoing flux 
in the asymptotic region. If to pass to non-stationary consideration of this question, then from the problem of 
possible singularity of the wave function at zero and sources we come to initial condition which should define 
further evolution of the Universe. In such a case, after defining the initial state (through set of parameters) it is 
possible to connect zero value of wave packet at a — (i. e. without singularity at such a point) with non-zero 
outgoing flux in the asymptotic region. In such direction different proposals have been made in frameworks of 
semiclassical models in order to describe inflation, to introduce time or to analyze dynamics of studied quantum 
system (for example, see [341135] ). 

In Ref. [15] an infinity potential wall at a = was introduced in order to consider the packets on the positive 
semiaxis of the scale factor a. A question can be appeared: how such an condition could be connected with 
possibility of the stationary wave functions to be non-zero at a = 0, according to logics above. In Appendix [Cl 
penetration of the packet is studied in fully quantum (non-semiclassical) consideration concerning the simplest 
barrier (which is enough to clarify this question). Here, we present exact analytical solutions for amplitudes. 
We find coefficients Tmir and Rmir describing penetration of the packet through the barrier from the internal 
region outside and its reflection, which are separated into the coefficients of penetrability Tb ar and reflection 

l 1 2 

i?bar of the barrier (in standard definition) and new coefficient \Ai nc \ describing oscillations of all multiple 
packets inside internal region: 

Tmir — \A- mc \ -Tbar, Rmir — \A- mc \ ■ i?bar- (37) 
This formulas seem to be fully quantum analog of the semiclassical formula of the decay width of a metastable 

l 1 2 

state introduced by Gurvitz and Kalbermann in Ref. |47) , where A lnc looks to be fully quantum analog of the 
normalization factor F. In Appendix [C] it is shown that consideration of the packet which penetrates through 
the barrier with its oscillations inside internal region provides complex stationary total wave function which is 
non-zero at a = also (while separate packets can be zero at such a point). So, we obtain full correspondence 
between non-stationary consideration of the packet penetrated through the barrier and described on the basis 
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of multiple internal reflections concerning boundaries, and the stationary formalism presented in this paper. 
Now difference between the penetrability calculated in Ref. jT5] and the penetrability obtained in this paper 
has became clear: in Ref. |15) the reflections of multiple packets are taken from boundary at a = into account, 
(in Ref. |15j Tmir is obtained, while I calculate Tb a r)- In order to pass from the rectangular barrier considered 
in Appendix C to arbitrary barrier shapes, formalism of Refs. [39,41,42 could be used, where we present exact 
solutions concerning barrier consisting from arbitrary finite number of rectangular steps of arbitrary sizes (which 
is supposed to approximate effectively the studied realistic barrier). 

5.5 The penetrability and reflection: the fully quantum approach versus semi- 
classical one 

Now we shall estimate by the method described above the coefficients of penetrability and reflection for the 
potential barrier <j3j) with parameters A = 36, B — 12 A, A = 0.01 at different values of the energy of radiation 
-Erad- We shall compare the found coefficient of penetrability with its value, which the semiclassical method 
gives. In the semiclassical approach we shall consider two following definitions of this coefficient: 

pWKB,(l) _ J_ pWKB,(2) 1 



penetrability n'2 ' penetrability / \ 2 ' f^S^ 

'20+l/(20) 2A 



where 



(BXt) 
t P 



exp J \V(a)-E\da. (39) 



int) 
tp 



One can estimate also duration of a formation of the Universe, using by definition (15) in Ref. [T3] : 

T = 2 dtp, int p " • (40) 

-^penetrability 

Results are presented in Tabl. [T] in Appendix 1. In calculations the coefficients of penetrability, reflection and 
mixing are defined by eqs. (|32|) . the fluxes by eqs. (|30|) . In the whole region of energies up to £ ra d = 2500 the 
calculations give coincidence of the first as a minimum 8 digits between -Pp^frabmty an< ^ ^pcncfrabiiity > therefore 
in the table the results for -Pp^frabmty are included only. From this table one can see that inside whole 
region of the energy the fully quantum approach gives value for the coefficient of penetrability enough close 
to its value obtained by the semiclassical approach, that differs essentially from results in the non-stationary 
approach [15j . This difference could be explained by difference in a choice of the boundary condition, which is 
used in construction of the stationary solution of the wave function. 

An advantage of the fully quantum method is possibility to calculate the coefficient of reflection. However, 
the calculations show that this coefficient inside the whole region of the energy used in Tabl. 1 equals to 1 
practically. But this coefficient differs visibly from 1 only at the energy of radiation enough close to the height 
of the barrier (see Tabl. 2 in Appendix), that is explained by essential decreasing of a role of the barrier in 
penetration of the wave through it. In order to estimate an accuracy of the found coefficients, I obtain property 
(|3"3"]) up to the first 11 digits (see Tabl. 2 in Appendix) inside whole region of the energy of radiation. The 
coefficient of mixing is less than 10 -19 (that is connected with computer error). So, there is no practical 
interference between the incident and reflected waves defined by such a way close to the internal turning point 
that points out their very accurate determination practically. Now it becomes clear that the approach proposed 
in Ref. [15] and the semiclassical methods do not give such an accuracy in the determination of the coefficients 
of penetrability and reflection. 



6 The penetrability in the FRW-model with the Chaplygin gas 

In order to connect the stage of universe with dust matter and its another accelerating stage, in Ref. |27j a new 
scenario with the Chaplygin gas was applied to cosmology. A quantum FRW-model with the Chaplygin gas has 
been constructed on the basis of equation of state instead of p (a) = p ra d(a)/3 (where p (a) is pressure) by the 
following (see also Refs. [251135] ): 

PCh = — — , (41) 
Pch 
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where A is positive constant and < a < 1. In particular, for the standard Chaplygin gas we have a = 1. 
Solution of equation of state (f4"Tj) gives the following dependence of density on the scale factor: 

/ B \ !/(!+<*) 

pch ( a )=^+-3 W J , (42) 

where B is a new constant of integration. This model through one phase a connects the stage of Universe 
where dust dominates and DeSitter stage. Note that for the first time the Chaplygin gas was introduced in 
aerodynamics [26) . 

Let us combine expression for density which includes previous forms of matter and the Chaplygin gas in 
addition. At limit a — > eq. (|42p is transformed into the pdust component plus the pa component. From such 
limit we find 

A = p A , B = p dust (43) 

and obtain the following generalized density: 



•2 , , 87r G I 2 ( Pdust \ . Prad I _ , ..s 



Now we have: 

3 V" ' a 3 ( 1+Q ) J ' a 2 (t) 

After quantization we obtain the Wheeler-De Witt equation 

l--^ + Vc h (a)\v(a)=E xad( p(a), E rild = (46) 

where 

«»W - (ire) "• J -2^" 4 (« + ^t) ■ (47) 

For the Universe of closed type (at k = 1) at 87rG = M p = 1 we have (see eqs. (6)-(7) in Ref. [55]): 

U C h(a) = 36a 2 -12a 4 (A + -^ T ) 1/(1+a) , i? rad - 12p rad . (48) 




Figure 7: Cosmological potentials with and without Chaplygin gas: Left panel is for potential V(a) = 36 a 2 — 
12 A a 4 with parameter A = 0.01 (turning point a tp = 17.320508 at zero energy E lad = 0), Right panel is for 
potential (|48|) with parameters A = 0.01, pdust = 30, a — 0.5 (minimum of the hole is -93.579 and its coordinate 
is 1.6262, maximum of the barrier is 177.99 and its coordinate is 5.6866). 

Let us expand the potential (|48[) close to arbitrary selected point a by powers q — a — a and restrict ourselves 
by linear item only: 

U C h (q) = V a + Vxq, (49) 
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For coefficients Vq and V± I find: 



Vo = Vch(a = a), 
Vi = dV Ch (a) 



da 



= 72 a 



U Ji{-4A-^}.(A. 



Pdust 
7 3(1+q) 



-a/(l+a) 



and eq. (|46[) obtains a form: 
After change of variable 



^2 ¥>(«) + (K) - S rad + Vi ?) = 0. 



eq. (|51|) is transformed into the following: 

d 2 



<K 



2^0 



(E^-Vo Vi \ 



After new change 



we have 



t _ E ia d — Vq Vi 

^ ~ |Vi| 2 /3 |Vi| C 



From eqs. ([5"2"j) and (|5"4")) we write: 



de 

E, 



v(0 + Cv(0 = o. 



rad 



Vo 



Vi 



3- 



(50) 

(51) 
(52) 

(53) 

(54) 
(55) 
(56) 



|Vi| 2 / 3 |Vi| 2 / 3 

Using such corrections after inclusion of the density component of the Chaplygin gas, I have calculated the 
wave function and on its basis the coefficients of penetrability, reflection and mixing by the formalism presented 
above. Now according to the logic of Sec. 14.11 I have defined the incident and reflected waves relatively new 
boundary which is located in the minimum of the hole in the internal region. Results are presented in Tabl. 3. 
One can see that penetrability is changed up to 100 times in dependence on the location of the boundary in 
such a coordinate or in the internal turning point (for the same fixed barrier shape and energy -E r ad)! This 
confirms that coordinate where incident and reflected waves are defined has essential influence on estimation 
of the coefficients of penetrability and reflection. This result shows real physical sense of the method proposed 
in the given paper. In the next Tabl. 4 I demonstrate fulfillment of the property (|33|) inside the whole energy 
region, which is calculated on the basis of the coefficients of penetrability, reflection and mixing obtained before. 
One can see that accuracy is the first 10-12 digits. Of course, any semiclassical calculations are not able to give 
such accuracy for penetrability for the studied cosmological barriers. 



7 Conclusions and perspectives 

In the paper the closed Friedmann-Robertson-Walker model with quantization in the presence of a positive 
cosmological constant and radiation is studied. I have solved it numerically and have determined the tunneling 
probability for the birth of an asymptotically deSitter, inflationary Universe as a function of the radiation 
energy. Note the following. 

1. A formalism for calculation of two linear independent partial solutions for the wave function of the Universe 
for the scale factor inside the region < a < 100 and the energy of radiation from zero up to the barrier 
height has been constructed. 

2. A fully quantum definition of the wave which propagates inside strong field and interact minimally with 
them, has been formulated for the first time, and approach for its stable determination has been con- 
structed. 

3. A new stationary approach for determination of the incident, reflected and transmitted waves relatively a 
barrier has been constructed, the tunneling boundary condition has been corrected. 
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4. A quantum stationary method of determination of coefficients of penetrability and reflection relatively 
the barrier with analysis of uniqueness of solution has been developed, where for the first time non-zero 
interference between the incident and reflected waves has been taken into account and for its estimation 
the coefficient of mixing has been introduced. 

5. A criterion of estimation of accuracy of the determination of these coefficients has been proposed on the 
basis of check of eq. ([33]) . 

In such a quantum approach the penetrability of the barrier for the studied quantum cosmological model with 
A = 36, B = 12 A parameters at A = 0.01 has been estimated with a comparison with results of other known 
methods. Note the following. 

• According to the calculations, inside whole region of energy of radiation the tunneling probability for 
the birth of an asymptotically deSitter, inflationary Universe is very close to its value, obtained in the 
semiclassical approach by eqs. (|3"5|) and but essentially differs on the results obtained before by the 
quantum non-stationary approach in Ref. [15] (see Tabl. 1 and 2 in Appendix). 

• The coefficient of reflection from the barrier in the internal region has been determined at first time. 
According to calculations, it is differed essentially on 1 at the energy of radiation close enough to the 
barrier height (see Tabl. 2 in Appendix). 

• The modulus of the coefficient of mixing is less 10 -19 for all energies, that points out that there is no 
interference between the found incident and reflected waves close to the internal turning point. 

• On the basis of the calculated coefficients I reconstruct a property (|3"3")l with accuracy of the first 11-18 
digits (see Tabl. 2 in Appendix) inside the whole studied region of the energy of radiation. 

One can assume that the method proposed can be easily generalized on the other cosmological models with the 
barriers with arbitrary complicated shapes. However, one can suppose that a visible change can be found in 
difference between the penetrabilities in the fully quantum and semiclassical approaches after inclusion of new 
components of density into the model (for example, after use of the component of Chaplygin gas in Refs. [26H29] ). 
Of course, it can be interesting to compare the estimations of the coefficients of penetrability and reflection 
with results which could be obtained on the basis of the Improved WKB approach on the basis of Refs. [30H32] . 



A Cosmological model in the Friedmann— Robertson— Walker metric 

A.l Dynamics of Universe in the Friedmann Robertson— Walker metric 

Let us consider a simple model of the homogeneous and isotropic Universe in Friedmann-Robertson-Walker 
(FRW) metric (see Ref. [T7], P- 438; also see Refs. [TSH5T]): 



where t and r, 9, <f> are time and space spherical coordinates, the signature of the metric is (—,+,+,+) as in 
Ref. |20) (see p. 4), a(t) is an unknown function of time and k is a constant, the value of which equals +1, 
or —1, with appropriate choice of units for r. Further, we shall use the following system of units: h = c = 1. 
For k = —1, the space is infinite (Universe of open type), and for k = +1 the space is finite (the Universe of 
closed type). For k = 1 one can describe the space as a sphere with radius a(t) embedded in a 4-dimensional 
Euclidian space. The function a(t) is referred to as the "radius of the Universe" and is called the cosmic scale 
factor. This function contains information of the dynamics of the expansion of the Universe, and therefore its 
determination is an actual task. 

One can find the function a(t) using the Einstein equations with taking into account of the cosmological 
constant A in this metric (we use the signs according to the chosen signature, as in Ref. [20] p. 8; the Greek 
symbols and v denote any of the four coordinates t, r, 9 and <f): 




(57) 



(58) 
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where is the Ricci tensor, R is the scalar curvature, T^„ is the energy-momentum tensor, and G is Newton's 
constant. From (|57p we find the Ricci tensor R^ v and the scalar curvature R: 

d ad a 2 bl Id 2 + ad + 2k 
Rtt = —3-, R rr = -r- + 2— — — = — — , 

a h n nr 1 — kr z , r „% 

h'r (59) 
R^ = R ee sin 2 6, Ree = ad r 2 + 2d 2 r 2 - h — + 1 = 2a 2 r 2 + ad 2 r 2 + 2£r 2 

r it d i „rr d , oe u , 00 p 6a 2 + 6qq + 6fc 
= .9 Rtt + 9 Rrr + 9 Ree + 9 R<P4, = ■ (60) 

The energy-momentum tensor has a form (see [20 , p. 8): — (p+p)U tl U^+pg^, where p andp are energy 
density and pressure. Here, one needs to use the normalized vector of 4- velocity U* = 1, U r = U 6 = = 0. 
Substituting the found components ([5^)) of the Ricci tensor i? M „, the scalar curvature ABU)) , the components of 
the energy-momentum tensor T^ v and including the component p r ad(a), describing the radiation in the initial 
stage (equation of state for radiation: p(a) — p ra( j(a)/3), into the Einstein's equation (|58[) at fj, = v = 0, we 
obtain the Friedmann equation with taking into account the cosmological constant (see p. 8 in Ref. [20] ; p. 3 in 
Ref. [21]; p. 2 in Ref. Q3]): 

where d is derivative a at time coordinate. From here, we write a general expression for the energy density: 

Prad 

a*(t) 



P(«)=PA + ^- (62) 



A. 2 Action, lagrangian and quantization 

We define the action as in Ref. [13] (see (1), p. 2): 



S = I V^9 [ - P ) dx 4 . (63) 



16ttG 

Substituting the scalar curvature (|60[) . then integrating item at d by parts with respect to variable t, we obtain 
the lagrangian (see Ref. [13] , (11), p. 4): 

c (a ' d) = £h ("" 2 + k ~ fl2 p(a) ) ■ (64) 

Considering the variables a and d as generalized coordinate and velocity respectively, we find a generalized 
momentum conjugate to a: 

d£(a,a) 3 . 

Pa = ^a— = -^G aa (65) 



1 [2-kG 2 2 3 fc 



and then hamiltonian: 

h(a,p a ) = pd-£(a,d) ^ - ^ |^^p 2 +a 2 ^^ -a 4 p(a)|. ( 66 ) 

The passage to the quantum description of the evolution of the Universe is obtained by the standard procedure 
of canonical quantization in the Dirac formalism for systems with constraints. In result, we obtain the Wheeler- 
De Witt (WDW) equation (see Ref. [13], (16)-(17), in p. 4, [UE1E2]), which after multiplication on factor and 
passage of the item at the component with radiation p ra d into right part transforms into the following form: 

{-5? + "«} *•>-*-*•>. v ^ = {ia) 2 ^-l^^ E —^- (67) 

This equation looks similar to the one-dimensional stationary Schrodinger equation on a semiaxis (of the variable 
a) at energy -E ra( j with potential V (a). For further analysis one can convenient to use the system of units where 
87r G = Af p = 1. Let us rewrite V (a) in a generalized form: 

V(a) = Aa 2 -Ba 4 . 



In particular, for the Universe of the closed type (k = 1) we obtain A = 36, B = 12 A (this potential coincides 
with [T5]). 



18 



A. 3 Potential close to the turning points: non-zero energy case 

At first, let us find the turning points at Pl i n and a tp , ut concerning the potential (|68[) at energy E,-^: 



A *BE iad A ABE r&d 

^^VM'f'f'T' atp ' out = V2B'V V ~ (69) 

Let us expand the potential V(a) Q68j) in powers q ou t = a — Ot p (where as dtp the point atp,in or (hp, out is used, 
close to which we find expansion), where (for small q) we restrict ourselves by the linear item only: 

V(q) = V + Vi, (70) 

where the coefficients Vq and V\ are: 

Vq = V(a = Otp, in) = V(ffl = Otp, out) = ^a? P _ - Ba t P = E *ad, 



- -Ws('-^)H>-^> 



finti , M / 4:BE rad \ f / 4B£ rad 



x 2B \ A 2 J \ V A 2 
Now eq. (j67[) transforms into a new form at variable g with potential V(q): 



dq 



-ZZ ¥>(<?) + Vi 9 y-fo) = 0. (72) 



B Calculations of the wave function of Universe 
B.l Method of beginning of the solution 

We shall be looking for the regular partial solution for the wave function close to arbitrary selected point a x . 
Let us write the wave function in the form: 

+oo H-oo 

<fi(a) = c 2 b n (a - a x ) n = c 2 b n a n , a = a - a x (73) 

and rewrite the potential through the variable a: 

V(a)^Co + C 1 d + C 2 d 2 + C 3 d 3 + C i a 4 , (74) 

where 

C = Aa 2 x -Ba x , 

Ci = 2a x {A- Ba 2 x ) -2Bal = 2Aa x -ABal, 

C 2 = A-Ba 2 x -ABa 2 x - B a 2 x = A-6Ba 2 x , (75) 
C 3 = -2Ba x — 2Ba x = — ABa x , 
d = -B. 

Substituting the wave function (|73[) . its second derivative and the potential (|74[) into Schrodingcr equation, we 
obtain recurrent relations for calculation of the unknown b n : 

_ (Cp - E) b _ (Cp -E)h+ C\ b _ (C -E)b 2 + d b 1 + C 2 b m) 
32 2 ' 3 6 ' 4 12 

, _ (Co - E) h + Cib 2 + C 2 bi + C 3 bp (77) 
5 20 ' 

(C - E) b n + d b n -t + C 2 b n - 2 + C 3 6„_ 3 + C 4 6„_ 4 , . . ,_ oN 

bn+2 = . w . at n > 4. 78 

(n + 1) (n + 2) 
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For given values for 60 and 61 using eqs. (|76|) — (fT5|) one can calculate all b n needed. At limit -EVad — ► and at 
a x = all found solutions for bi are transformed into the corresponding solutions (40), early obtained in |16j at 
-Erad = 0. Using c 2 = 1, from eqs. ([75)1 we find: 

6 = ip(a x ), bi = ip'(a x ). (79) 

So, on the basis of the coefficients bo and b\ one can obtain values for the wave function and its derivative at 
point a x . Implying two different boundary conditions through bo and b\ (in such a way that they locate the 
first node for the wave function at different places) , we obtain two linearly independent partial solutions tpi (a) 
and if2 (a) for the wave function. Using the internal turning point at P) i n as the starting point, by such a way we 
calculate the first partial solution which increases in the barrier region (we select: 60 = 0.1, b\ = 1), and using 
the external turning point dt P) out as the starting point, by such a way we calculate the second partial solution 
which decreases in the barrier region (we select: b = 1, &i = —0.1). Such a choice provides effectively a linear 
independence between two partial solutions. 



B.2 Method of continuation of the solution 

Let us rewrite eq. (15?) in such a forrrjj 

ip"(a) = f(a)<p(a). (80) 

Let {a„} be a set of equidistant points a n — ao + nh. Denoting values of the wave function ip (a) at points a„ 
as tp n , we have constructed own algorithm of the ninth order for determining <^ n +i and tp' n on the previously 
known tp n and tp n -i- 



3ii+5oi . 3oi 3io -Sooffn . n ,,9^ 

ip n+ l = tp n -i \-tp n hO ft , 

3oi - 3n 3oi ~ 3n 

/ 2 . 3io - 3oo . ^ /, 

<Pn = fn-1 H fn h O (h ), 

3oi - 3n 3oi - 3n 



where 

3oo 



(81) 



2+h\ fn + ^ (f::+f n ) + ±h e (/i 4) +4(/;) 2 + 7 / n /:+/3)+ 
+ I h& (ti 6) + 16 fn ; " 4) + 26 /; ^ + 15 ^") 2 + 22 ; " ; « + 28 fn ^") 2 + ; " 

301 = l^ + l^/^ + G^/^ (82) 

310 = § h 3 f n + 1 h 5 (f^ + 4 /„/;) + 1 h 7 (/ w + 11 /„ fW + 15 /x + 9 f n /;) , 

311 = 2 h + 1 h 3 f n + 1 ^ 5 (3 + fi) + j [ h 7 (h / w + 13 /„ + 10 (/;) 2 + /^) . 

A local error of these formulas at point a n equals to: 

5n = ^h w f^. (83) 



C Tunneling of the packet through radial rectangular barrier 

Let us consider a problem of quantum tunneling of the packet through the barrier used in cosmological model. 
We shall study such a process consequently by steps of its propagation relatively to each boundary of the barrier, 
using developed formalism of multiple internal reflections presented in Refs. [?.? J3"6H4"0] (see also Refs. [43H46] ). 
In order to form idea of multiple internal reflections of the packets in description of quantum tunneling on the 
positive semiaxis of the scale factor a, we shall use the simplest potential V(a): V(a) = —Vq for < a < R\ 
(internal region I), V{a) = V\ for R x < a < R 2 (region II of the barrier) and V(a) = for a > R 2 (external 
region III). For simplicity, we start from consideration of the case when total energy of system E is higher then 
the barrier height V\ : E > V\ . 

2 in obtaining this algorithm a logics from 1231 was used 
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In the first step we consider the packet in the region I propagating to the right, which is incident on the first 
boundary of the barrier at R\: 

ip^ c {a,t) = J g(E - E)e lkia - lEt ' h dE at < a < R 1} (84) 

-Emin 

where k± = \J E + Vq , E is the energy. The weight amplitude g(E — E) can be used in standard form of gaussian 
and satisfies to normalization J \g(E — E)\ 2 dE — 1, value E is an average energy. This packet transforms into 
two new packets: the first packet transmitted through this boundary and propagating further in the region II, 
and the second one reflected from the boundary and propagating back in the region I: 

ipll\a,t) = J g{E ~ E)a^ e la2a - lEt / h dE at i?i < a < R 2 , 

^ „ (85) 

^ll(a,t) = / g{E- E)A { x> e- lkia - lEt l n dE at < a < R u 

Emin 

where k 2 = y/E — V\ . We find new unknown coefficients and , using requirements of continuity of the 
total wave function ip(a,t) (which is summation of all packets) and its derivative at R±: 

a (l) = _^*^ e «*i-fa)*i, A$ =h^l e 2ikiR K (86) 
fci + k 2 ki + k 2 

In the second step we consider further propagation of the packet ip[l\a,t), which is incident on the second 
boundary at R 2 . It transforms into two new packets: the first packet transmitted through this boundary and 
propagating in the region III, and the second one reflected from this boundary and propagating back in the 
region II. We define these packets in the form 

V> (2) (M) = + f g(E-E)tpW(a)e- lEt / h dE, (87) 

Emin 

where as the stationary parts we use: 

Vinc(a) = a (1) e 4fc2a , for R 1 < a < R 2 , 

^(a) = A^e ika , fora>i? 2 , (88) 
Vref( a ) = /3 (1) e^ fc2a , for Ri < a < R 2 , 

where k — y/E. Imposing condition of continuity on the total wave function and its derivative at R 2 , we obtain 
two new equations, from which we find new unknowns coefficients Aj? and : 

A™ = T+ ■ a« , T+ = e* R - , 

/3« = R+ • Q W, R+ = ^-Jle 2ik - R \ 
2 2 k 2 + k 

We have introduced two new coefficients T 2 , R 2 , which logically connect the transmitted and reflected ampli- 
tudes Aj, ' and with the incident amplitude in this step. Here, we shall use bottom index for denotation 
of number of the considered boundary, upper (top) sign "+" or "— " for positive (to the right) or negative (to 
the left) direction of the incident wave, correspondingly. So, we can write T± — and Rf = A^ also 

,,(2) 



R 

In the third step we consider further propagation of the reflected packet ip^tl in the region II. Incidenting 

(3) 

on the first boundary, it transforms into new packet iptr > transmitted through this boundary and propagating 

(3) 

in the region I, and into new packet Vvef j reflected from boundary and propagating back in the region II. We 
define the new packets by eq. (|87p (with upper index 3), where as the stationary parts we use: 

= ^cf(«)> f or Ri < a < Ba, 
tp£\a) = A ( ^e- ikia , farO<a<-Ri, (90) 
(p®j,{a) = a^e ik2a , for R 1 < a < R 2 . 
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From continuity conditions for the total wave function and its derivative at Ri , we find the unknowns coefficients 
A% } and aW: 



a (2) =R ~ .0(1), R~ = , _, 



i (fei fej) fii 



ki + k 2 (qd 
fci + A; 2 

In the forth step we need to consider further propagation of the reflected packet ip^~l in the region I in 
the 1-st step. It is incident on the first boundary at a = transforming into new packet propagated to the 
right. At such a point we can include different considerations of origin of possible sources at a = 0, possible full 
propagation (like in spherically symmetric problems of quantum decay in nuclear physics which is 3-dimensional 
and we have no additional boundaries at a = 0) or, in contrary, full reflection used in different fully quantum 
approaches (like introduction of an infinite potential wall at a — in Ref. [H]). In order to produce ability to 
work with different such considerations, we write: 

<^(a) = ^(a), <p^(a,k 1 )=R^-^tl(a,-h)=Aitle ik ^, for < a < R lt (92) 

where 

4tl=Ro-^l (93) 

Supposing the full propagation through this boundary (without any possible reflections), we obtain Rq = — 1. 
If we liked to use the condition of the infinite potential wall at a — 0, than we should also have Rq = — 1. 

Analyzing further reflections and transmission of the packets concerning the boundaries, we conclude that 
any of following steps is similar to one of 4 considered above. From analysis of these steps recurrent relations 
are found for calculation of new unknown amplitudes A^\ A^' and for arbitrary step n, 

summations of these amplitudes are calculated. These series can be calculated easier, using coefficients 7^ and 
Rf. Analyzing all possible "paths" of the propagations of all possible packets inside the barrier and internal 
well, we obtain: 

+oo , +0O 1 

E At' = i + RfRo+ Rt Ro -RfRo + -. = 1+ E (Ri V 



= 1 m = l ^ ' ' 1 — Ri Rq 

e4° = (e 4f)-K n+T+.R+R-,- r+ + ...} = (e41)-^ 

71=1 V n=l ' y ' x n=l ' 



A^ — Fty + Fty ■ R R^ + R± ■ Rq R± ■ R R± + ... — 



71 = 1 



71=1 



where 



K = ^ n m t, • [ i + e {Ri Ry) m ) = Ri + 



m=l ' 1 -rt-2 /ggN 

+oo _ mN T.+ T + \ ' 



^ = T+T+.(l+ + f (R+R-) m )=-^^ 



m—l 

The resultant expressions for the incident, transmitted and reflected packets concerning the barrier are 
written in form of eq. (|87|) . where the following stationary wave functions should be used: 

^inc(a) = e lkia , for < a < Ri, 

+00 

ptr(a) = E ^e ika , for a>R 2 , 

71=0 \^°) 

-\-oo 

(p re{ (a) = J2 A n R e~ ikia , for0<a<i?i. 

71=0 

At finishing, we determine the full amplitudes 

+00 . . +00 f s +00 /i™ +00 

AT=E4 n) , 4=E4 B) ! «=E« (n) = ^, = E = " • i? 2 + (97) 

n— 1 n— 1 n— 1 2 n— 1 



22 



and coefficients T and R describing penetration of the packet from the internal region outside and its reflection 
from the barrier ^ ^ 

TmIR = T" \At\ = \A- mC \ ' Jbar, ^bar = T~ 1^1 , . , 

«i «i (98) 

RMIR = \Ar\ — \A inc \ ■ Rbar, Rbar = \R>l\ > 

where Tb ar and i?bar are coefficients of penetrability and reflection of the barrier (in standard definition), and 

l 1 2 

A- lnc is coefficient determining oscillations of the packet inside the internal region (this is fully quantum analog 
of the normalization factor F introduced in Ref. [47] for semiclassical description of nuclear decay). 

Series J2 A< £ } C , J2 a t \ E a r\ E« (n) and Y,P {n) obtained using the approach of the multiple internal 
reflections, exactly coincide with the corresponding coefficients A- mc , At, Ar, a and f3 calculated by standard 
stationary method (where the continuity conditions of the stationary total wave function and its derivative are 
used at each boundaries, and the wave function is not equal to zero at a = 0). We test property: 

\A T \ 2 + \A R \ 2 = 1 or T M ir + Rmir = 1, (99) 

ki 

which is fulfilled and confirms that the method MIR gives us proper solution for the wave function. If energy is 
less then the height of the barrier, then for description of penetration of the wave through such a barrier with 
its tunneling it needs to use the following change [37ll39ll40] : 



k 2 ^i£, Z = y/E-Vi. (100) 

Using it, all found above solutions are applied for the problem with tunneling through the barrier. For the 
barrier consisting from two rectangular steps of arbitrary heights and widths we have already obtained exact 
coincidence between amplitudes for the wave function obtained by method of MIR and the corresponding 
amplitudes obtained by standard approach of quantum mechanics. Increasing of number of the rectangular 
steps in the barriers keeps such a coincidence and fulfillment of the property (|99|) (see some generalizations in 
Refs. [39,41,42 ). So, we have obtained full coincidence between all amplitudes calculated by method MIR and 
by standard approach of quantum mechanics, and that is way we generalize the method MIR for description 
of penetration of the wave through barrier consisting from arbitrary number of rectangular steps of arbitrary 
sizes. 



D Calculations of the penetrability 

The penetrabilities of the barrier with A — 36, B — 12 A parameters at A = 0.01 for the closed Friedmann- 
Robertson-Walker model with quantization in the presence of a positive cosmological constant and radiation 
calculated by the fully quantum approach and by the semiclassical approach are presented in Tabl. 1. One 
can see that inside the whole region of the energy of radiation .E^ad up to 2690 the penetrability calculated by 
the fully quantum method is very close to its value, obtained by the semiclassical approach by eqs. Q38p and 
(|39"1) . Also one can find an essential difference between these calculations and results obtained in Ref. [TS] by 
the non-stationary quantum approach. In the next Tabl. 2 the coefficients of the penetrability, reflection and 
mixing calculated in the fully quantum method are presented for the energy of radiation E ra< \ close to the height 
of the barrier. One can see that summation of all such values for coefficients allows to reconstruct the property 
(f3"3"|) with accuracy of the first 11-18 digits. 
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Energy 
E la d 
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7.66 


15.53 


1800.0 


6.4255 x 
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1.2437 x 
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2.4783 x 
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1.2803 x 
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7.96 


15.38 
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1.1189 x 
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2.1580 x 
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1.4776 x 
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7.6619 x 
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8.26 


15.22 


2000.0 


1.3288 x 
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2.5653 x 
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1.2914 x 
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7.3137 x 10" 
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13.37 
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1.6663 x 10- 
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3.0428 x 10" 
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11.07 


13.31 
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6.9240 x 10" 
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1.2606 x 10" 
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1.7678 x 10+ 14 


11.14 


13.25 
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2.8842 x 10" 


-11 


5.2116 x 10" 
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4.3050 x 10+ 12 


11.21 
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1.2002 x 10 


-9 


2.1495 x 10" 
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1.8825 x 10+ 10 


1.0511 x 10+ 11 


11.29 


13.12 


2650.0 


4.9881 x 10 


-8 


8.8401 x 10 


-9 


4.5642 x 10+ 8 


2.5754 x 10+ 9 


11.38 


13.05 


2660.0 


2.0738 x 10 


-6 


3.6263 x 10 


-7 


1.1068 x 10+ 7 


6.3303 x 10+ 7 


11.47 


12.97 


2670.0 


8.7110 x 10 


-5 


1.4836 x 10 


-5 


2.6596 x 10+ 5 


1.5615 x 10+ 6 


11.58 


12.87 


2680.0 


3.6953 x 10 


-3 


6.0519 x 10 


-4 


6.3369 x 10+ 3 


3.8693 x 10+ 4 


11.70 


12.76 


2690.0 


1.5521 x 10 


-1 


2.4634 x 10 


-2 


1.5293 x 10+ 2 


9.3602 x 10+ 2 


11.86 


12.61 



Table 1: The penetrability -Ppcnctrabmty of the barrier and the duration r of the formation of the Universe defined by 
eq. (|40[) in the fully quantum and semiclassical approaches 
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Energy 


Fully quantum method 


Turning points 


-Erad 


Penetrability 


Reflection 


Mixing 


Summation 
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"■tp, out 


2690.0 


0.15521440329121 


0.84478559670782 


1.47 x 


1(T 


19 


0.99999999999904 


11.86 
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2691.0 
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1.59 x 


10" 


19 


0.99999999999871 
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12.59 
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7.61 x 


io- 


20 


0.99999999999831 


11.90 


12.57 


2693.0 


0.42120001498898 
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1.41 x 


10- 


19 


0.99999999999785 


11.93 


12.55 


2694.0 


0.55773509442073 
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1.44 x 


io- 


19 
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11.95 
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7.82 x 


10" 


20 


0.99999999999701 


11.98 
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7.80 x 
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3.45 x 


10" 
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0.99768097743782 
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2.10 x 
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12.07 
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1.43 x 


10" 


20 


0.99999999999901 


12.12 


12.36 



Table 2: The coefficients of the penetrability, reflection and mixing calculated by the fully quantum method and their 
summation 



Energy 


Penetrability P pen etrability 




Time 


T 


Full QM method 1 


Full QM method 2 


Method WKB 


Full QM method 1 


Method WKB 
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-4 


3.4471 x 10+ 5 


1.2530 x 10+ 4 


170.0 


3.4894 x 10 


-3 


2.6078 x 10' 


-1 


4.2919 x 10" 


-2 


3.0460 x 10+ 3 


2.4820 x 10+ 2 



Table 3: The penetrability Ppenetrability of the barrier and duration r of the formation of the Universe defined by eq. (I40[) 
in the FRW-model with the Chaplygin gas obtained in the fully quantum and semiclassical approaches (minimum of the 
hole is -93.579 and its coordinate is 1.6262, maximum of the barrier is 177.99 and its coordinate is 5.6866): the fully QM 
method 1 is calculations by the fully quantum approach for the boundary located in the coordinate of the minimum of 
the internal hole (i. e. coordinate is 1.6262), the fully QM method 2 is calculations by the fully quantum approach for 
the boundary located in the internal turning point dtp, in (coordinates of the turning points are in Tabl. 4) 
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Table 4: The coefficients of the penetrability, reflection and mixing calculated by the fully quantum method and test on 
their summation for the FRW-model with the Chaplygin gas density component (the fully quantum approach 1 is used 
at the internal boundary located in the coordinate of the minimum of the internal hole) 
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